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Abstract 

The responses of perturbations added into the optically thick, advection-dominated 
accretion disk (ADAD), what we call the slim disk (SD), are investigated through 
numerical simulations. Although it is proposed that the SD is thermally stable, I find 
that a perturbation added into the disk is not rapidly damped and moves through 
the disk in its free-fall time. After the perturbation moves, the global structure 
of the disk does not vary very much. These facts may account for the substantial 
variability of the X-ray luminosities of stellar super-luminal jet sources (SLJSs) and 
Narrow-Line Seyfert Is (NLSls). 
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1 Introduction 

Recent X-ray observations report that not only stellar black hole candidates 
(SBHCs) in their low states (=faint sources), but also NLSls and stellar SLJSs 
(=bright sources) exhibit X-ray fluctuations (variability). The fluctuations of 
bright sources are made in the optically thick ADAD, what we call SD (=bright 
disk). However, the time evolution of the SD has not been investigated so 
far. The numerical simulation by Manmoto et al. (1996) is well known as 
a time-evolution calculation of the optically thin ADAD (=faint disk). The 
disturbance added into the optically thin ADAD falls into the central star. 
The disk luminosity increases when the disturbance falls, and the light curve 
of this process is in good agreement with the X-ray shot configuration of the 
SBHC Cyg X-l in its low state. 

I investigate how the luminosities of disturbed SDs vary. I add a similar dis- 
turbance as Manmoto et al. (1996) into the SD, and as a result, a similar light 
curve to that of the optically thin ADAD is obtained. 

The basic equations and numerical procedures are described in §2. The resul- 
tant time evolution and discussion will be presented in §3. 
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Table 1 

Parameter sets of Models A-C. 



Model 


ro 


k 


A 


20 


-0.3 


B 


20 


-0.1 


C 


50 


-0.3 



2 Basic Equations 



I calculate the evolution of a one-dimensional axisymmetric disk. The basic 
equations are the same as those of Manmoto et al. (1996), and are those of 
mass conservation, momentum conservation, angular momentum conservation, 
and energy flow: 
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where E(= / pdz) is the surface density, W(= J pdz) is the vertically inte- 
grated pressure, and e is the internal energy of the accreting gas. Q(— v^/r) 
and f2 K [= (GM/r) 1 / 2 / (r— rs)] are the angular frequency of the gas flow and the 
Keplerian angular frequency in the pseudo-Newtonian potential (Paczyhski & 
Witta 1980), respectively, where M is the mass of the central black hole and 
rs is the Schwarzschild radius. I set the viscosity parameter to be a = 0.1. 

To evaluate the radiative cooling rate, F, I consider black body radiation: 

F ~ 37R/2 + V/3' (5) 



where a is the radiation constant, T is the gas temperature, and r R is the 
Rosseland optical depth. 

I obtain the steady state solution of equations 1-4 and perform the time- 
evolution calculation using the steady state solution as the initial state. As for 
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the outer boundary condition, I set all quantities to be those of the standard 
disk by Shakura & Sunyaev (1973). The inner boundary is set at r in = 2.7rs, 
where a free boundary is adopted. I add a mass of MAt into the disk through 
the outer boundary at every time step (At). I set M — 100L Edd /c 2 where L Edd 
is the Eddington luminosity. 



I add a perturbation (disturbance) to the initial state of the disk as 

bp 
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Throughout these calculations, I assign the wavelength of the perturbation, 
A = 20rs- I set two parameters, the radius of the center of the perturbation, 
r , and the ratio between 5p and p at r = r , k. These parameters are var- 
ied, with three calculations being performed. The parameter sets of the three 
calculations are listed in Table 1. 



3 Results and Summary 



Fig. 1 presents the light curves of Models A-C. All models show shot-like light 
curves. 

Fig. 2 plots how the surface density, S, varies after the perturbation is added 
for Model A. The dotted and solid lines represent the initial state and the 
later evolution, respectively. The perturbation does not change its configura- 
tion very much, and it does not rapidly decay. The time for the perturbation 
to propagate corresponds to the free-fall time. After the perturbation has 
damped, the disk structure is not globally changed. 

The response of SDs to local disturbances has been examined by one-dimensional 
numerical simulations. It is generally believed that SDs are thermally stable. 
I, however, find that disturbances added into the accretion flow do not damp 
rapidly and decay with roughly the free-fall time. After the disturbance has 
damped, the global disk structure of the disk is not greatly modified. This can 
account for the persistent X-ray emission with substantial variations observed 
in NLSls and SLJSs. When a perturbation is made in the SD, it decays and 
exhibits one X-ray shot. Since the structure of the SD does not globally vary 
much after a perturbation propagates, the next perturbation produced in the 
disk can also exhibit an X-ray shot. Repeating such processes continuously 
can make the substantial variability seen in the X-ray luminosities of SLJSs 
and NLSls. 
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Fig. 1. Light curves of Models A-C. The solid lines represent the light curves, and 
the dotted line represents the luminosity in the steady state. 
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Fig. 2. The time evolution of the disk surface density, S. The dotted and solid 
lines represent the initial state and the later evolution, respectively. The solid line 
labeled 1 means the initial configuration of the disturbance. Elapsed times for the 
labels 2-7 are t = 397, 785, 1352, 2894, 5997 and 14288r s /c, respectively. Note 
that both axes represent the relative values. 
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